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Abstract 

> 

Consideration of the analytic properties of pion-induced baryon self-energies 

^^ . leads to new functional forms for the extrapolation of light baryon masses. 

These functional forms reproduce the leading non-analytic behavior of chiral 
perturbation theory, the correct non-analytic behavior at the Ntt threshold 
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ON 

ON ' and the appropriate heavy-quark limit. They involve only three unknown 



parameters, which may be obtained by fitting to lattice data. Recent dy- 
namical fermion results from CP-PACS and UKQCD are extrapolated using 
Q-i, these new functional forms. We also use these functions to probe the limit 

r- 1 of applicability of chiral perturbation theory to the extrapolation of lattice 

QCD results. 
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I. INTRODUCTION 

In the last year there has been tremendous progress in the computation of baryon masses 
within lattice QCD. Improved quark |l[ and gluon || actions, together with increasing 
computer speed means that one already has results for N, A and vector meson masses for 
full QCD with two flavors of dynamical quarks. Although the results are mainly in the 
regime where the pion mass (m^) is above 500 MeV, there has been some exploration as low 
as 300 - 400 MeV on a 3.0 fm lattice by CP-PACS [|. 

In spite of these impressive developments it is still necessary to extrapolate the calcu- 
lated results to the physical pion mass (yU = 140 MeV) in order to make a comparison with 
experimental data. In doing so one necessarily encounters some non-linearity in the quark 
mass (or m^), including the non-analytic behavior associated with dynamical chiral sym- 
metry breaking. Indeed, the recent CP-PACS study did report the first behavior of this 
kind in baryon systems. 

As the computational resources necessary to include three light flavors with realistic 
masses will not be available for many years, it is vital to develop a sound understanding of 
how to extrapolate to the physical pion mass. We recently investigated this problem for the 
case of the nucleon magnetic moments || . 

The cloudy bag model (CBM) U is an extension of the MIT bag model incorporating 
chiral symmetry. It therefore generates the same leading non-analytic (LNA) behavior as 
chiral perturbation theory (xPT). This model was recently generalized to allow for variable 
quark and pion masses in order to explore the likely mass dependence of the magnetic 
moment ||. This work led to several important results : 



• 



a series expansion of /%,( n ) in powers of m n is not a useful approximation for m^ larger 
than the physical mass, 

on the other hand, the behavior of the model, after adjustments to fit the lattice data 
at large m^ was well determined by the simple Pade approximant: 

Eq. (HD not only builds in the Dirac moment at moderately large m\ but has the correct 
LNA behavior of chiral perturbation theory 

jj, = /Jo — am n , 

with a a model independent constant, as m% — > 0. 

fixing a at the value given by chiral perturbation theory and adjusting /io and /3 
to fit the lattice data yielded values of /x p and //„ of 2.85 ± 0.22 hn and —1.96 ± 
0.16/ijv, respectively, at the physical pion mass. These are significantly closer to the 
experimental values than the usual linear extrapolations in m q . 
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FIG. 1. One-loop pion induced self-energy of the nucleon and the delta. 

Clearly it is vital to extend the lattice calculations of baryon magnetic moments to lower 
values of m n than the 600 MeV used in the study just outlined. It is also important to 
include dynamical quarks. Nevertheless, the apparent success of the extrapolation procedure 
suggested by the CBM study gives us strong encouragement to investigate the same approach 
for baryon masses. 

Accordingly, in this paper we study the variation of the N and A masses with m n (or 
equivalently m q ). Section [II| is devoted to considerations of the low- lying singularities and 
pion-induced cuts in the complex plane of the nucleon and A spectral representation. The 
analytic properties of the derived phenomenological form are consistent with both chiral 
perturbation theory and the expected behavior at large m q . This phenomenological form 
is eventually fitted to recent two-flavor, full QCD measurements made by CP-PACS || 
and UKQCD f/j. However, to gain some insight into the parameters and behavior of the 



functional form we examine the N and A masses as described in the CBM in section III 



In section |V| we apply the analytic form to the lattice data. Section |V] is reserved for a 



summary of our findings. 



II. ANALYTICITY 



By now it is well established that chiral symmetry is dynamically broken in QCD and 
that the pion is almost a Goldstone boson. As a result it is strongly coupled to baryons 
and therefore plays a significant role in the N and A self energies. In the limit where the 
baryons are heavy, the pion induced self-energies of the N and A, to one loop, are given by 
the processes shown in Fig. |T|. Note that we have restricted the intermediate baryon states 
to those most strongly coupled, namely the N and A states. 

The analytic expression for the pion cloud correction to the masses of the N and A are 
of the form || 

5M N = a NN + a NA , (2) 

where 

„ 3 2 / j. k u NN (k) 

CTNN = ~ 1fi 2f2 9A / «« 2 n\ ' ( S > 

iSTr 2 /^ Jo w 2 (k) 

3 32 2 j k u NA (k) 

aNA -~16^fl25 gA J Q d w(k)(AM + w(k))' {A) 



and 

5 M A = o"aa + ctatv , (5) 

where 

0"AA = (T NN , (6) 

— ^ ^ „2 f°° A u k m JVa(^) f 7 \ 

aAN ~ T&^fl2B 9A J d w(k)(AM - w(k)) ■ [V 

We note that AM = M& — M^ , qa = 1-26 is the axial charge of the nucleon, w(k) = 
^Jk 2 + m 2 is the pion energy and v,NN(k), UNA(k), . ■ . are the NNn, NAir, . . . form factors 
associated with the emission of a pion of three momentum k. We have used SU(6) symmetry 
to relate the four coupling constants to the NNn coupling, which, in turn, has been related 
to gA/2fn by chiral symmetry. The form factors arise naturally in any chiral quark model 
because of the finite size of the baryonic source of the pion field - which suppresses the 
emission probability at high virtual pion momentum. As a result, the self-energy integrals 
are not divergent. 

The leading non-analytic contribution (LNAC) of these self-energy diagrams is associated 
with the infrared behavior of the corresponding integrals - i.e., the behavior as k — > 0. As a 
consequence, the leading non-analytic behavior should not depend on the details of the high 
momentum cut-off, or the form factors. In particular, it should be sufficient for studying the 
LNAC to evaluate the self-energy integrals using a simple sharp cut-off, u{k) = 9(A — k). 
In Section |TTTj we shall compare the results with those calculated using a phenomenological, 
dipole form factor and show that this is in fact an effective simplification. 

Using a ^-function for the form factors, the NNtt and AA7r integrals (c.f. Figs. |l|(a) and 
|l](d), respectively), which are equal, are easily evaluated in the heavy baryon approximation 
used here: 
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The integral corresponding to the process shown in Fig.|](b), with a #- function form factor, 
may be analytically evaluated. For m n > AM 

9 2 a , , 2 a 2^/2 f f Jml + A 2 + AM + A 

"™ = -55f% (12( ™" - AM)> { aICtm ( - ^ral-^P 

AM + m n \ 
il-AM 2 ) 



— arctan 
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+3AM(3ml - 2AM 2 ) In ( V^±E±A 

\ m w 



3 V^ 2 + A 2 AM A + 6 AM 2 A - 6m 2 A + 2A d ) , (9) 



while for m w < AM we find 
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(10) 



Similar results are easily obtained for the process shown in Fig.[l](c). For m v > AM, the 
analytic form is 
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(12) 



The self-energies involving transitions of N — > A or A — > N are characterized by the branch 
point at m^ = AM. 



A. Chiral Limit 

Chiral perturbation theory is concerned with the behavior of quantities such as the 
baryon self-energies as m q —* 0. For the expressions derived above, this corresponds to 



taking the limit m^ — > 0. The leading non- analytic (LNA) terms are those which corre- 
spond to the lowest order non-analytic functions of m q - i.e., odd powers or logarithms of 
mn. By expanding the expressions given above, we find that the LNA contribution to the 
nucleon/Delta mass (Eq.®) is given by 

M L N ?£ = -^ff>l, (13) 

in agreement with a well known result of x?T || . A careful expansion of the An contribution 
to the nucleon self energy, Eq.(0), yields the LNA term: 

3o 2 32 3 
a NA (m„ A) ~ ^f^SAM^ ln( ^ } (14) 



as m v — > which is again as expected from x?T [E03- F° r the Nir contribution to the 



self-energy of the A, the LNA term in the chiral limit of Eq. ([11]) yields: 

a AN (m„A) ~ ~^fA^M<Hrn,) ■ (15) 

Of course, our concern with respect to lattice QCD is not so much the behavior as 
m n — > 0, but the extrapolation from high pion masses to the physical pion mass. In this 
context the branch point at m 2 = AM 2 is at least as important as the LNA near m n = 0. 



We shall return to this point later. We note that Banerjee and Milana jlT] found the same 
non-analytic behavior as m„- — > AM that we find. However they were not concerned with 
finding a form that could be used at large pion masses - i.e. one that is consistent with 
heavy quark effective theory. 

B. Heavy Quark Limit 

Heavy quark effective theory suggests that as m^ — > oo the quarks become static and 
hadron masses become proportional to the quark mass. This has been rather well explored 
in the context of successful nonrelativistic quark models of charmonium and bottomium 
|12]| . In this spirit, corrections are expected to be of order l/m q where m q is the heavy 
quark mass. Thus we would expect the pion induced self energy to vanish as l/m q as the 
pion mass increases. The presence of a fixed cut-off A acts to suppress the pion induced self 
energy for increasing pion masses, as evidenced by the m 2 in the denominators of Eqs. (^|), 
(|) and (0). While some m 2 dependence in A is expected, this is a second-order effect and 
does not alter the qualitative features. By expanding the arctan (A/m„-) term in Eq. @ for 
small A/run, we find 

3ffj A 5 Ln fA 7 \ 
aNN --TQ^l5^ + °WJ ' (16) 

which vanishes for m n — > oo. Indeed, in the large m^ (heavy quark) limit, both Eqs. (|9p 
and (|Tl|) tend to zero as 1/m 2 ,. 
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C. Analytic form 

We now have the chiral and heavy quark limits for each of the four integrals in Fig.|I|. 
These expressions, which contain a single parameter, A, are correct in the chiral limit - i.e., 
they reproduce the first two non-analytic terms of xPT. They also have the correct behavior 
in the limit of large pion mass, namely they vanish like l/m^. The latter feature would be 
destroyed if we were to retain only the LNA pieces of the self-energies as they would diverge 
at large m n faster than m q . Rather than simplifying our expressions to just the LNA terms, 
we therefore retain the complete expressions, as they contain important physics that would 
be lost by making a simplification. 

We note that keeping the entire form is not in contradiction with ;\TT, as we have 
already shown that the leading non-analytic structure of ;\TT is contained in this form. 
However as one proceeds to larger quark (pion) masses, differences between the full forms 
and the expressions in the chiral limit will become apparent. For example, the branch 
point at m^ = AM 2 , which is an essential non-analytic component of the m^-dependence 
of the self-energy and which should dominate in the region m n ~ AM, is also satisfactorily 
incorporated in Eqs.(|J) and flTT|). Yet the LNA chiral terms given in section |il A| know 
nothing of this branch point and are clearly inappropriate in the region near and beyond 
m 2 = AM 2 . 

As a result of these considerations, we propose to use the analytic expressions for the 
self-energy integrals corresponding to a sharp cut-off in order to incorporate the correct 
LNA structure in a simple three-parameter description of the m^- dependence of the N and 
A masses. In the heavy quark limit hadron masses become proportional to the quark mass. 
Moreover, as we shall see in the next section, the MIT bag model leads to a linear dependence 
of the mass of a baryon on the current quark mass far below the scale at which one would 
expect the heavy quark limit to apply. This is a simple consequence of relativistic quantum 
mechanics for a scalar confining field. On the other hand, lattice calculations indicate that 
the scale at which the pion mass exhibits a linear dependence on m q is much larger than that 
for baryons.Q In fact, over the range of masses of interest to us, explicit lattice calculations 
show that m 2 is proportional to m q . Hence we can simulate a linear dependence of the 
baryon masses on the quark mass, m q , in this region, by adding a term involving m 2 . The 
functional form for the mass of the nucleon suggested by this analysis is then: 

M N = a N + (3 N ml + a NN (m n , A) + cx/va^, A) , (17) 

while that for the A is: 

M A = «a + PauiI + ctaa(^tt, A) + a an (m*, A) . (18) 

The mass in the chiral limit is given by 



x One doesn't expect such linear behavior to appear for quark masses lighter than the charm quark 
mass where the pseudoscalar mass is 3.0 GeV. Even at this scale the quarks are still somewhat 
relativistic. 



M^ = a N + a NN (0,A)+(T NA (0,A), (19) 

where the meson cloud effects are explicitly contained in o/vat(0, A) + ct/va(0, A). The mass 
of the A in the chiral limit is calculated in an analogous way. We know that flTT| ) and flTSp 
have the correct behavior in the chiral limit. Individually, they also have the correct heavy 
quark behavior .[] Between the chiral and heavy-quark limits there are no general guidelines, 
so in the next section we shall compare our functional form to the Cloudy Bag Model, a 
successful phenomenological approach incorporating chiral symmetry and the correct heavy 
quark limit. 



III. BARYON MASSES WITHIN THE CBM 

As a guide to the quark mass dependence of the N and A masses we consider the Cloudy 
Bag Model (CBM) [J6"Ul3"f. This is a minimal extension of the MIT bag model such that chiral 
symmetry is restored, which has proven quite successful in a number of phenomenological 
studies of baryon properties and meson-baryon scattering |||15|-|17|] . Within the CBM, a 



baryon is viewed as a superposition of a bare quark core and bag plus meson states. The 
linearized CBM Lagrangian with pseudovector pion-quark coupling (to order 1/ f n ) is 



C=[q{i<fdr-m q )q-B]6v-~ 

+ 1 -(d^) 2 - \mW + %j-Jtfl*rq ■ d,n , (20) 

where B is the bag constant, f n is the n decay constant, 9 V is a step function (unity inside 
the bag volume and vanishing outside) and 5s is a surface delta function. In a lowest order 
perturbative treatment of the pion field, the quark wave function is not effected by the pion 
field and is simply given by the MIT bag solution fl9|^|2Tf . 

In principle the tcNN form factor can be directly calculated within the model. It dies 
off at large momentum transfer because of the finite size of the baryon source. Rather 
than using this calculated form factor, which is model dependent, we have chosen to use a 
common phenomenological form, namely a simple dipole: 

where k is the magnitude of the loop (3-)momentum, /i is the physical pion mass (139.6 
MeV), and A D is a regulation parameter. 



2 With regard to the difference, M/\ — Mn, HQET suggests that this difference should vanish as 
m n — > oo. This is only guaranteed by Eqs. ( |T7D and (O) (through Eq. (|T6|)) if the entire mass 
difference arises from the pion self-energy. While one could enforce this condition through the 
introduction of additional parameters and a more complicated analytic structure for the higher- 
order terms of Eqs. ([17]) and (|i~8|), we prefer to focus on the regime of m 2 from 1 GeV 2 to the 
chiral limit. As we shall see, Eqs. ( [171 ) and (^) are quite adequate for this purpose. 
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FIG. 2. The pion mass dependence of the N and A baryons generated in the CBM using a 
dipole form factor with A^ = 1 GeV. Fits of fll7|) and ( |l8|) to the CBM results are illustrated by 
the curves. 



In the standard CBM treatment, where the pion is treated as an elementary field, the 
current quark mass, m q , is not directly linked to m n . Most observables are not sensitive to 
this parameter, as long as it is in the range of typical current quark masses. For our present 
purpose it is vital to relate the m q inside the bag with m n . Current lattice simulations 
indicate that m^ is approximately proportional to m q over a wide range of quark masses 
0. Hence, in order to model the lattice results, we scale the mass of the quark confined 
[m^l ' jj) m q , with m q being the current quark mass corresponding to 



in the bag as m q 

the physical pion mass \x 



rn 



(0) 



is treated as an input parameter to be tuned to the lattice 



results, but in our magnetic moment study it turned out to lie in the range 6 to 7 MeV, 
which is very reasonable. 

The parameters of the CBM are obtained as follows. The bag constant B and the 
phenomenological parameter z$ are fixed by the physical nucleon mass and the stability 
condition,[] dM N /dR = 0, for a given choice of R and m q ' . For each subsequent value 
of the pion mass or the quark mass considered, u> and R are determined simultaneously 
from the linear boundary condition [|T9|--|2T|| and the stability condition. In this work we 
have calculated the mass of the iV and A baryons as a function of squared pion mass (as 
illustrated in Fig. [2]). The A calculation is similar to that for the N, however the value of B 
is fixed to be the same as that used for the nucleon, while z is adjusted to fit the observed 
mass difference, taking into account the pionic contribution to this quantity, at the physical 



3 Note that while zq,B and the ttNN form factor may all depend on m q , this dependence is ex- 
pected to be a smaller effect and we ignore such variations in order to avoid an excess of parameters. 
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FIG. 3. Pion induced self-energy corrections for a 1 GeV dipole form factor. The LNA term of 
^PT tracks the NNir contribution up to m^ ~ 0.2 GeV, beyond which the internal structure of 
the nucleon becomes important. 
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As expected on quite general grounds (and discussed in Section [IT]), as the pion mass 
increases the mass of the baryon does indeed become linear in m^. In addition, from the 
curvature at low pion mass, we see that the non-analytic structure is important in the region 
m n below 400 MeV. 

We now fit our functional forms for the baryon masses, Eqs. (|i~7f) and fll8|) , to the CBM 
data. We note that the CBM data is generated using a phenomenologically motivated, 
dipole form factor, whereas the functional form used in the fit involves a 9 cut-off. In order 
to simulate the fitting procedure for lattice data, our fit to the CBM results involves only 
pion masses above the physical branch point at Ma = M^, followed by an extrapolation to 
lower pion mass. 

It can be seen from Fig. |2| that our extrapolation to the physical pion mass is in good 
agreement with the CBM calculations: at the physical pion mass the extrapolated A" mass 
is within 0.8% of the experimental value to which the CBM was fitted, while the A is within 
0.3% of the experimental value. We present the parameters of our fit in Table |. The value 
for the sharp cut-off (A) is 0.44(2) GeV, compared to A^ =1 GeV for the dipole form factor. 

It was noted in Section |T| that the constant a in our functional form is not the mass of 
the baryon in the chiral limit, but rather this is given by M N = aN + aNN(0,-h) + 0"iVA(O, A) 
- with an analogous expression for the A. We find that the extrapolated A" and A masses 
in the chiral (SU(2)-flavor) limit are (M^\m^) = (905,1210) MeV, compared with the 
CBM values (898,1197) MeV. 

The mass dependence of the pion induced self-energies, a^, for the 1 GeV dipole form 
factor, is displayed in Fig. H. The choice of a 1 GeV dipole corresponds to the observed 



axial form factor of the nucleon |22] , which is probably our best phenomenological guide to 



10 



the pion-nucleon form factor p3| . We note that Ujvjv tends to zero smoothly as m n grows 
and it is only below m 2 ~ 0.3 GeV 2 that there is any rapid variation. That this behaviour 
cannot be well described by a polynomial expansion is illustrated by the dotted curve in 
Fig. |[ There we expanded ctjvat about m n = as a simple polynomial, a + /3m 2 + jm^, 
with 7 fixed at the value required by chiral symmetry. Clearly the expansion fails badly for 
m n beyond 300-400 MeV. 

The behavior of the A% contribution to the self-energy of the A is especially interesting. 
In particular, the effect of the branch point at m n = AM is seen in the curvature at 
m 2 ~ 0.1 GeV 2 . For comparison, we note that while there is also a branch point in the 
nucleon self-energy at the same point - see Eq. (^) - the coefficient of (m 2 — AM 2 ) 2 vanishes 
at this point. As a consequence there is little or no curvature visible in the latter quantity 
at the same point. The correct description of this curvature is clearly very important if one 
wishes to obtain the AN mass difference at the physical pion mass. The fact that, as shown 
in Fig. 0, our simple three parameter phenomenological fitting function can reproduce N 
and A masses within the CBM, including this curvature, suggests that this should also 
provide a reliable form for extrapolating lattice data into the region of small pion mass. 

Fig. [| illustrates the degree of residual model dependence in our use of (17) and (18). 
There the variation of the nucleon self-energy, <Jnn, calculated with a 1 GeV dipole form 
factor (solid curve) is fit using the form a + /?m 2 +<jNN(m 7T , A) (dash curve, with a = —0.12 
GeV, j3 = 0.39 GeV -1 and A = 0.57 GeV). Note that the deviations are at the level of a 
few MeV. For the A the self-energy, a an, is again calculated using a 1 GeV dipole form 
factor and fit with our standard fitting function, a + /3ra 2 + OA^m^, A). The quality of the 
fit (with a = -0.062GeV, (3 = 0.024 GeV" 1 and A = 0.53 GeV) is not as good as for the 
nucleon case. Nevertheless, the difference between the two curves at the physical pion mass 
(vertical dotted line) is only about 20 MeV. At the present stage of lattice calculations this 
seems to be an acceptable level of form factor dependence for such a subtle extrapolation. 

IV. LATTICE DATA ANALYSIS 

We consider two independent lattice simulations of the N and A masses, both of which 
use improved actions to study baryon masses in full QCD with two light flavors. The 
CP-PACS |3| lattice data is generated on a plaquette plus rectangle gauge action with im- 
provement coefficients based on an approximate block-spin renormalization group analysis. 
The (9(a)-improved Sheikholeslami and Wohlert clover action is used with a mean- field im- 
proved estimate of the clover coefficient cgw=l.QA - 1.69. This estimate is likely to lie low 



Baryon 


a 

(GeV) 



(GeV- 1 ) 


A 

(GeV) 


M B 
(GeV) 


Error 


N 

A 


1.09 
1.37 


0.739 
0.725 


0.455 

0.419 


0.948 
1.236 


0.8% 
0.3% 



TABLE I. Parameters for fitting (O) and (18) to the CBM data. Here we have taken 
Ro = 1.0 fm and m q = 6.0 MeV. The Error column denotes the relative diff< 
experimental values which were used as a constraint in generating the CBM data. 
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FIG. 4. Comparison between the nucleon and A self-energies, ctnn and ctaa^ calculated us- 
ing a dipole form factor (solid and long-dash dot curves, respectively) and fits using the form 
a + (3m^ + aij(m n , A), based on a sharp cut-off in the momentum of the virtual pion (dash and 
short-dash dotted curves respectively). 



relative to a nonperturbative determination |L4] and may leave residual O(a) errors. 

Ideally one would like to work with lattice data in which the infinite-volume continuum 
limit is taken prior to the chiral limit. Until such data is available, we select results from 
their 12 3 x 32 and 16 3 x 32 simulations at /3 = 1.9. Lattice spacings range from 0.25 fm to 
0.19 fm and provide physical volumes 2.7 fm to 3.5 fm on a side. While the volumes are 
large enough to avoid significant finite volume errors, the coarse lattice spacings necessitate 
the use of improved actions. Systematic uncertainties the order of 10% are not unexpected. 

The UKQCD [0] group uses a standard plaquette action with the (9(a)-improved Sheik- 
holeslami and Wohlert action. At a /3 of 5.2, UKQCD uses csw — 1-76, which is lower 
than the current non-perturbative value [14j] of 2.017, again leaving some residual 0(a) er- 
rors. Lattice spacings are necessarily smaller, ranging from 0.13 to 0.21 fm. We select their 
12 3 x 24 data set providing better statistical errors than their largest volume simulation. 
Physical volumes are 1.6 to 2.6 fm on a side suggesting that finite volume errors may be an 
issue on the smallest physical volume where the dynamical quark mass is lightest. 

In full QCD, the renormalized lattice spacing is a function of both the bare coupling and 
the bare quark mass. In order to determine the lattice spacing, the UKQCD collaboration 
calculates the force between two static quarks at a distance r$ |24|, while CP-PACS considers 
the string tension directly. While the two approaches yield similar results in the quenched 
approximation, string breaking in full QCD may introduce some systematic error in the 
extraction of the string tension at large distances. In fact we find that the two data sets are 
consistent, provided one allows the parameters introducing the physical scale to float within 
systematic errors of 10%. A thorough investigation of these systematic errors lies outside 
the scope of this investigation. Instead we simply rescale the UKQCD and CP-PACS data 
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Scaling 




TV 








A 






CP-PACS 


UKQCD 


a 





M N 


i< 


a 





M A 


M™ 






(GeV) 


(GeV" 1 ) 


(GeV) 


(GeV) 


(GeV) 


(GeV" 1 ) 


(GeV) 


(GeV) 


0% 


0% 


1.10 


0.778 


0.954 


0.910 


1.29 


0.680 


1.150 


1.125 


+5% 


-5% 


1.15 


0.736 


1.003 


0.961 


1.36 


0.602 


1.227 


1.203 


0% 


-10% 


1.10 


0.767 


0.957 


0.914 


1.31 


0.624 


1.169 


1.145 


+10% 


0% 


1.20 


0.707 


1.050 


1.008 


1.42 


0.581 


1.285 


1.262 



TABLE II. Parameters for fits of ( |l7| ) and (|l~8|) to lattice data. Here we fix A (A^r = 0.455 and 
Aa = 0.419) and vary a and 0. The mass of the baryon at the physical pion mass is Mn (M/\) 
and the mass in the chiral limit is ML (M* ). The scaling columns represent adjustments to the 
scale parameters providing physical dimensions to the lattice data. 

sets in combining them into a single, consistent data set. 

We begin by considering the functional form suggested in Section [TJ with the cut-off A 
fixed to the value determined by fitting the CBM calculations. The resulting fits to the 
baryon masses are shown in Fig. ||] for the unshifted lattice data and Fig. [| where each data 
set is adjusted by 5% to provide consistency. The extrapolations are indicated by the solid 
(dashed) curve for N (A). The resulting fit parameters and masses^ are listed in Table |J. 

In examining fits in which the cutoff is allowed to vary as a fit parameter, we found it 
instructive to also study the dependence of the fit on the number of points included. This 
dependence is shown for the N in Fig. |7] and for the A in Fig. |^. In particular, we 
compare fits including the lowest lattice point (at around 0.1 GeV 2 ) and then excluding it. 
When we fix the value of A the fits are stable and insensitive to the lowest point. They tend 
to lie slightly above the lowest data point. However, given the caution expressed by the 
CP-PACS collaboration for the lowest point, we view these fits as reasonably successful. In 
contrast, when the value of A is treated as a fitting parameter, it is sensitive to the inclusion 
of the lowest point. Hence, to perform model independent fits, it is essential to have lattice 
simulations at light quark masses approaching m 2 ~ 0.1 GeV 2 . An analysis of the current 
data suggests A = 0.661 GeV and provides a nucleon mass 130 MeV lower than the CBM- 
constrained fit. Tables |TTT| and [TV] summarize parameters and physical baryon masses for N 
and A respectively. 

It is common practice in the lattice community to use a polynomial expansion for the 
mass dependence of hadron masses. Motivated by %PT the lowest odd power of m n allowed 



is ml: 



M N = a + (3ml + l m l (22) 



The results of such fits are shown in Figs, pi and |10| for N and A respectively. The corre- 



sponding parameters are reported in table [V[ As can be seen in table [V], the coefficient of 



4 The errors bars for the extrapolated baryon masses at the physical pion mass displayed in the 
figures are naive estimates only. We are unable to perform a complete analysis without the lattice 
results on a configuration by configuration basis. 
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FIG. 5. Baryon masses calculated by UKQCD (open points) and CP-PACS (solid points), as 
a function of m^. The solid (dashed) curve illustrates a fit to the combined data sets for N (A). 
The left-most data points are our extrapolated values of the baryon masses at the physical pion 
mass. 
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FIG. 6. UKQCD and CP-PACS baryon masses with 5% adjustments in the scale parameters 
to improve the agreement between the two data sets. (The key is as described in Fig. [51) 
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FIG. 7. UKQCD and CP-PACS nucleon masses with scale parameters adjusted by 5%. The 
data is as described in Fig. || The dashed lines represent fits without the point at 0.1 GeV 2 . The 
solid lines include this point. The top pair of lines are fits with A fixed at 0.455 GeV, a value 
preferred on the basis of our CBM analysis. The bottom pair have A as a fit parameter. 
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FIG. 8. UKQCD and CP-PACS A-baryon masses with scale parameters adjusted by 5%. The 
data is as described in Fig. || The dashed lines represent fits without the point at 0.1 GeV 2 . The 
solid lines include this point. The top pair of lines are fits with A fixed at 0.419 GeV, a value 
preferred on the basis of our CBM analysis. The bottom pair have A as a fit parameter. 
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The two parameter fit corresponds to using Eq. 22 with 7 set equal to the value known from xPT. 



The three parameter fit corresponds to letting 7 vary as an unconstrained fit parameter. The solid 
line is the fit for the functional form of ([17]) , fit (d) of Table ID . 
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is the fit for the functional form of (18), fit (d) of Table IV. 
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the m^ term, which is the leading non-analytic term in the quark mass, disagrees with the 
coefficient known from ;\PT by almost an order of magnitude. This clearly indicates the 
failings of such a simple fitting procedure. We recommend that future fitting and extrap- 
olation procedures should be based on Eqs. fll7|) and (|i~8l) , which are consistent with x?T 
and the heavy quark limit. 

V. SUMMARY 

In the quest to connect lattice measurements with the physical regime, we have explored 
the quark mass dependence of the N and A baryon masses using arguments based on 
analyticity and heavy quark limits. In the region where m n is larger than 500 MeV, the 
lattice data can be reasonably well described by the simple form, a + /3mJ, which is linear 
in the quark mass. The additional curvature associated with chiral corrections only appears 
below this region. This can be understood quite naturally within chiral quark models, like 
the cloudy bag, which lead to a cut-off on high momentum virtual pions, thus suppressing 
the self-energy diagrams quite effectively as m^ increases. The pionic self-energy diagrams 
which we consider are unique in that only these diagrams give rise to the leading non- 
analytic behaviour which yields a rapid variation of baryon masses in the chiral limit. Loops 
involving heavier mesons or baryons cannot give rise to such a rapid variation. 

Based on these considerations, we have determined a method to access quark masses 
beyond the regime of chiral perturbation theory. This method reproduces the leading non- 
analytic behavior of x?T and accounts for the internal structure for the baryon under 
investigation. We find that the predictions of the CBM, and two flavor, dynamical fermion 
lattice QCD results, are succinctly described by the formulae of equations (|T7|) and fll8D with 
terms defined in (||) through (|T2|). We believe that equations (||) - (0) are the simplest one 
can write down which involve a single parameter, yet incorporate the essential constraints 
of chiral symmetry and the heavy quark limit. 

Firm conclusions concerning agreement between the extrapolated lattice results and ex- 
periment cannot be made until the systematic errors can be reduced below the current level 
of 10% and accurate measurements are made at m^ ~ 300 MeV or lower. The significance 
of nonlinear behavior in extrapolating nucleon and A masses as a function of m% to the 
chiral regime has been evaluated. We find that the leading non-analytic term of the chi- 



Fit a (3 
(GeV) (GeW 1 ) 

(a) L76 0.386 

(b) 1.15 0.727 

(c) 1.42 0.564 

(d) 1.15 0.736 

TABLE III. Parameters for the fits shown in Fig. [7]. Parameter sets (a) and (b) are obtained 
by excluding the lowest data point from the fit, while (c) and (d) include it. Parameter sets (a) 
and (c) are fits with 3 parameters, sets (b) and (d) are fits with A fixed to the phenomenologically 
preferred value. 
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A 


M N 


(GeV) 


(GeV) 


0.789 


0.763 


0.455 


1.010 


0.661 


0.870 


0.455 


1.003 



Fit a (3 A M A 
(GeV) (GeV 1 ) (GeV) (GeV) 

(a) L64 0.414 0.683 1.042 

(b) 1.37 0.587 0.419 1.240 

(c) 1.54 0.475 0.616 1.095 

(d) 1.36 0.602 0.419 1.230 

TABLE IV. Parameters for the fits shown in Fig. |8|. Parameter sets (a) and (b) are obtained 
by excluding the lowest data point from the fit, while (c) and (d) include it. Parameter sets (a) 
and (c) are fits with 3 parameters, sets (b) and (d) are fits with A fixed to the phenomenologically 
preferred value. 

ral expansion dominates from the chiral limit up to the branch point at m^ = AM. The 
curvature around m n = AM, neglected in previous extrapolations of the lattice data, leads 
to shifts in the extrapolated masses of the same order as the departure of lattice estimates 
from experimental measurements. 
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Fit 


a 





7 or A 


M N 


a 


(3 


7 or A 


M A 




(GeV) 


(GeV- 1 ) 


(GeV^ 2 ) or (GeV) 


(GeV) 


(GeV) 


(GeV- 1 ) 


(GeV" 2 ) or (GeV) 


(GeV) 


(a) 


-0.128 


7.38 


-5.60 


-0.001 


0.182 


7.09 


-5.60 


0.304 


(b) 


0.912 


1.69 


-0.761 


0.943 


1.18 


1.45 


-0.703 


1.202 


(c) 


1.15 


0.736 


0.455 


1.003 


1.37 


0.602 


0.419 


1.227 



TABLE V. Parameter sets for the fits shown in Figs. ^ and [K]. Set (a) is for the 2 parameter 
fit of ( P2[) with 7 from xPT, (b) for the 3 parameter fit of (P2|), and (c) for the preferred functional 
form. 
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